Strange Heat Flux in (An)Harmonic Networks 
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Abstract 

We study the heat transport in systems of coupled oscillators driven out of equilibrium 
by Gaussian heat baths. We illustrate with a few examples that such systems can exhibit 
"strange" transport phenomena. In particular, circulation of heat flux may appear in the 
steady state of a system of three oscillators only. This indicates that the direction of the heat 
fluxes can in general not be "guessed" from the temperatures of the heat baths. Although we 
primarily consider harmonic couplings between the oscillators, we explain why this strange 
behavior persists under weak anharmonic perturbations. 
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1 Networks of oscillators 

In this note, we consider steady states of (an)harmonic oscillators driven by heat reservoirs at 
different temperatures. We show, by simple examples, that "anything is possible" for such physi- 
cal systems: in particular, it is basically impossible to guess in which direction energy flows. We 
will first describe the harmonic case and then argue why the results extend to mildly anharmonic 
problems. 



The setup is that of n masses, all equal to 1, connected by a set of harmonic "springs," at 
most n{n — l)/2 of them. For the sake of simplicity, the position and velocity of each mass are 
chosen to be one-dimensional. The potential is a function V(qi, q 2 , . . . , q n ), which is given as a 
positive definite quadratic form i(q, V^q). The (Gaussian) heat baths interact with some (at least 
2) of the n masses. Each mass is either attached to its own heat bath at temperature T, L > 0, 
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with friction T.;, or is attached to no heat bath. In this case we will say that 1^ = and leave Tj 
undefined. The stochastic differential equations describing such a system are for i — 1, . . . , n: 



dpi = -(Vq)i dt - TiPi dt + \J2TiTi duji{t) , 
dqi = pidt , 

where the uJi(t) are independent Wiener processes. It will be convenient to write the problem in 
matrix form. Let x = (p, q) denote the state of the 2n masses. The invariant measure of the 
problem (if it exists), is (up to normalization) of the form exp(— |(x, Q _1 x)) with the 2n x 2n 
matrix Q being the solution to the Lyapunov equation 

QA* + AQ = -B , (1) 

where 

f-T -V\ n (2YT 

A =(i o J ■ B ={ 

Here, T and T are the diagonal matrices whose elements are Tj and Tj. We denote by H = {i : 
Ti ^ 0} the indices of the masses in direct contact with a heat bath. The following condition 
assures uniqueness of the invariant measure, and can be easily derived from [1] : 

Lemma 1 Consider the space S spanned by the vectors {V k ei, i G TC, k = 0, . . . , n}, where et 
denotes the i th unit vector of R™. If S = R n , then (1) has a unique solution. Moreover, this 
solution is positive definite. 

When the condition of Lemma 1 is not satisfied, a change of coordinates shows that at least 
one mode is neither coupled to a heat bath nor to the rest of the system. The simplest exam- 
ple where this happens is shown in Fig. 1 (see [2, 3]). The masses 1 and 2 are coupled to heat 
baths, while the masses 3 and 4 are only coupled to the masses 1 and 2. All the springs have the 
same coupling constant. Writing the equations of motion, one easily checks that the variables 
q = q 3 — g 4 and p = p 3 — p 4 evolve as an isolated harmonic oscillator. 




Figure 1 : Non-unique steady state 



We henceforth assume that the assumptions of Lemma 1 hold (this can be easily verified for 
the examples given in the sequel). Therefore, the steady state exists, is unique, and we denote 
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by (/) the average value of an observable f(x) in this state. For convenience, we shall write the 
matrix Q as four n x n blocks 




where X and Y are positive symmetric matrices. As a consequence of (1), R* = —R. Averages 
of quadratic observables are given by the elements of the matrix Q, namely (piPj) = Xij, 
(liQj) = Y ij and (PiQj) = 



2 Heat fluxes 

We briefly recall a common definition of a heat flux between two points of the system [4]. In 
general, the evolution of an observable / is given by the equation / = Cf, where C is the 
Fokker-Planck operator, in our case 

C = p ■ V 9 - q ■ V V p - p ■ TV P + V p ■ YTV P . 

By definition we have (Cf) = 0. The energy in the spring connecting points i and j is Uij = 
— hVij(qi — qj) 2 , where V^- < when the coupling is attractive. In order to obtain the heat flux 
between these two points, we interpret the equation (CUij) = as a conservation equation for 
the energy in the spring, and identify the terms in this equation as energy fluxes. We denote the 
average value of the flux from i to j by (f)^, whose expression turns out to be 

= Vij (pj(qi - qj)) = Vij (pjqi) , 

since = Ra = by antisymmetry of the matrix R. For a point i connected to a bath, the 

heat flux entering the system through that point, denoted by i5 is obtained similarly, leading to 

4>i = Tj (Tj — {p\)) . (3) 

Because of energy conservation, the total heat flux at every point has average zero in the steady 
state. In the sequel, we only consider average quantities and by flux we always mean average 
flux in the steady state. 



Very few results are available concerning the direction of the heat fluxes in the system. The 
main one is the positivity of the global entropy production, namely 

This (strict) inequality has been proved in [4] for an anharmonic chain between two baths at 
different temperatures. Under the conditions of Lemma 1, one can easily show that it remains 
valid for harmonic networks [2]. Because the matrix Q is positive definite, we can also conclude: 
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Lemma 2 The point(s) attached to the hottest bath(s) cannot absorb heat from the other baths. 
The point(s) attached to the coldest bath(s) cannot inject heat in the system. 

Proof. If all the temperatures are the same, say one easily checks that the steady state is 
Gibbsian, that is 

1 



« = *U v") ■ (4) 

In this equilibrium state all the fluxes vanish and the lemma is trivially verified. Consider next 
a system S with at least two different temperatures, and denote by Q the solution to the corre- 
sponding Eq.(l). Let T max be the temperature of the hottest heat bath(s) of S and > T max be 
an arbitrary higher temperature. We define a system S' as a copy of S but whose temperature 
matrix X" is given by 



r 



Ti > for % G H , 
otherwise . 



Let Q' be the solution to Eq.(l) for S'. We note that when all the parameters but T are fixed in 
Eq.(l), the solution Q(T) is linear in T. Therefore Q + Q' is a Gibbsian matrix (4) with $ = 0, 
in particular 

x« + x , a = e, 

where we have used the block notation (2) for Q and Q' . Since both matrices are positive- 
definite, we have Xu and X' u > 0, therefore X u < for any > T max , and finally Xu < T max . 
We consider next the flux <^ entering the system through a "hot" point i for which Tj = T max . 
Because of Eq.(3) we have 

</>i = Ti (T max - (pi)) = T t (T max - Xu) > . 

The corresponding inequality for the cold point(s) is obtained by an equivalent construction. This 
concludes the proof of the Lemma 2. 

The two results we have mentioned give some information on how the system of oscillators 
exchanges heat with the baths. We are now interested in knowing how the flux propagates within 
the system of oscillators. The main observation of this note is that "everything" is possible, 
basically through superposition of elementary solutions. Indeed, by the linearity of Q(T), each 
heat bath can be considered as an independent flux source, and the total flux at any point is simply 
obtained by adding the contribution of all the baths. This is how the four examples below can be 
found. 

Example 1 A linear chain 



Consider a linear chain composed of four equal masses, each of which is coupled to a heat bath. 
In the setup of Fig. 2, the heat flux is going against the (local) temperature gradient between the 
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Ti = 100 



To 



T, = 20 



Ta 



(1)juuuuuuuuuuuu(^^^ 
Figure 2: The flux between 2 and 3 goes against the temperature gradient 



l,Vl2 = V 23 = V34 



13.5,02^3 = 1.1,03- 



2.8) 



masses 2 and 3. Instead of defining the local temperature with the heat bath, we can also use 
the local kinetic energy (pi). However, we find (pi) = 14 < 18 = (p§), thus the "backward 
flux" persists. This first example can be easily understood: as noted in the proof of Lemma 2, 
the matrix Q solving (1) is a linear function of T and so are the fluxes. Thus, we can decompose 
our system as the sum of two similar chains, one with temperatures (100, 0, 0, 2) and the other 
with (0, 2, 20, 0). The total flux in the middle spring still goes to the right, since temperature Ti 
pushes much more energy into the chain than T 3 does. 

Example 2 Circulation of heat 




Figure 3: A circulation of energy remains in the steady state 

(Ti = T 2 = 1, V12 = -10, V13 = V 23 = -20, 0i_ 3 = 02-^3 = 0.290, 0i^2 = 0.008) 



In this second example, the heat injected in the system by the hot bath has two possible "chan- 
nels" to reach the cold bath. What is surprising is the appearance of a "backward flux" in one of 
them which is not due to excess temperature as in Example 1. As a result of this, a circulation 
of heat remains in the steady state, as shown in Fig. 3. This example shows that energy fluxes 
between heat baths, as we understand them in this note, are not similar to electrical currents be- 
tween potentials. Indeed, should the arrows of Fig. 3 represent electrical currents, the potentials 
Ui at points i — 1, 2 and 3 should satisfy U\ > U2 > U3 > U\. In other words, Fig. 3 contra- 
dicts a "Kirchoff's Law" on current loops. Such an example can also be constructed when the 
"triangle" is in the center of a chain connecting two heat baths. 

Example 3 Three connected heat baths with different coupling constants 
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T 2 = l 



T 3 = l 



Figure 4: Circulation in a "fully thermalized" triangle 

(Ti = 1, V 12 = Viz = -45, V13 = -30, 0i^2 = 0.57, <fc_ 3 = 0.35, <fe^l = 0.03) 



The example in Fig. 4 shows that the circulation of Example 2 can also be produced when all 
three masses are in contact with a bath, even if two baths have the same temperature. 

Example 4 A "heat pump" 



In this last example, we construct a system that mimics a thermodynamic heat pump. Figure 5 
shows two chains of three oscillators coupled through their middle point. The ends of the upper 
chain are connected to the hottest (Ti) and the coldest (T 2 ) bath, while the ends of the lower one 
are connected to intermediate temperatures (T 3 and T 4 ). Here again, the heat in the lower chain 
flows against the temperature gradient. The interesting point in this example is that no energy is 
flowing (in average) between the two chains (05^6 = 0). It is as if the upper chain was acting 
on the lower one through fluctuations only. By slightly varying the temperatures, one can even 
obtain 05_>6 < 0. This last variant is quite different from the thermodynamics of Carnot cycles, 
since the subsystem in which heat is pumped releases energy into the pump. We remark that 




4>T 4 = 6.8 



iyT 2 = 3.6 



Figure 5: The lower part of the system pumps heat from cold to hot 

(r* = i, Vis = v 56 = v 36 = 40, y 25 = v 46 = 20, <^ 5 = 2.4, ^ 6 ^ 3 = 0.2, <fe^ 6 = o> 
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Lemma 2 prevents us from building a pump between two baths if one of them is an extremal 
temperature of the system. 

The situations described in these four examples are not a "far from equilibrium" behavior, 
because there is no such thing in harmonic systems. Indeed, because of the linearity of Q(T), the 
fluxes in the system keep the same sign when (all) temperatures are rescaled. Moreover, fluxes 
are unchanged when all the temperatures are shifted by a constant. Therefore, everything we have 
shown can happen at arbitrary temperatures, and with arbitrarily small temperature differences. 
We also remark that because of the loops of heat flux, there is no possible definition of local 
temperature that would prevent the flux from going against the temperature gradient. 



3 Weakly anharmonic systems 

It is well-known that the heat transport in harmonic systems does not reproduce the usual macro- 
scopic laws, in particular Fourier's law does not hold [5]. Justification is commonly seen in the 
fact that the modes of a harmonic system are extended, causing the heat to be transported bal- 
listically rather than diffusively. Note however that phenomena described in this note continue 
to hold in a weak anharmonic limit. Indeed, consider a slight perturbation of the coupling, for 
instance 

Ve(q) = l(q,Vq)+ £ -J2c lJ (q i -q J ) i . 

The existence of a unique steady state for certain systems with such a potential is proved in [6]; 
since every point in our examples is reached in a simple way from a heat bath, the results of [6] 
generalize to this case. The corresponding (smooth) invariant measure p £ is not Gaussian but 
still decays rapidly at large energies. If p £ as a function of e is sufficiently regular around e = 0, 
fluxes are continuous functions of this parameter. Then, for every example we have shown, one 
can find a small enough e so that fluxes of the perturbed system have the same direction as those 
in the unperturbed system. Although a written proof of the regularity of p £ does not seem to be 
available yet, this result is believed to be true [7, 8]. A key point is the following: as explained 
in [9], the system in the steady state spends most of its time below a certain energy level, with 
only rare excursions to high energies. With sufficiently small e, one can make sure that the an- 
harmonicity is irrelevant in arbitrary long parts of the dynamics. 
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